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Abstract. We discuss differences in shapes, expansion 
velocities and fragmentation times of structures created 
by an energy deposition from a single Gamma-Ray Burst 
(GRB) and an OB association to the ISM. After the initial 
inflation, supershells produced by GRBs are almost static 
or slowly expanding, contrary to more rapidly expanding 
supershells created by OB associations. We discuss the po- 
sition of the energy source relative to the symmetry plane 
of the galaxy: observed arc-like structures can be the most 
dense part of structures formed by an expansion from a 
source above or below the galactic plane. Arcs may also 
form if the expansion takes place inside a giant HI cloud. 
We try to reproduce the size, the age, and the average 
distance between OB associations in the Sextant region 
at the edge of LMC 4. 

Stars: formation ISM: bubbles ISM: su- 
pernova remnants Galaxies: individual: LMC 
Galaxies: star clusters Gamma rays: bursts 



1. Introduction 

It is generally assumed that the HI shells and supershells 
observed in galaxies may be formed by the action of mas- 
sive stars in OB associations or they may be created by 
impacts of high velocity clouds (HVC) into the galac- 
tic HI plane. Recently, another possibility has been pro- 
posed: they may result from an energy input connected 
with Gamma-Ray Burst (GRB) events. This may follow 
coalescence of two neutron stars in a binary system lead- 
ing to the formation of a black hole. Such events, last- 
ing a few milliseconds, release up to 10 54 erg of orbital 
kinetic energy of the progenitor. The energy may be re- 
leased into a collimated jet illuminating only a fraction 
fb of its sky. Such beaming would decrease the output 
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energy by a factor /& and increase the frequency of occur- 
rence by a factor fZ~ . Another way to release several times 
10 52 erg of energy is the collapse of a single super-massive 
(> 1OOM ) star, the hypernova (Paczynski 1998). X-ray 
afterglows, discovered by the satellite BeppoSAX, lasting 
several hours, result from the interaction between rela- 
tivistic particles and the ambient medium around a black 
hole. The discovery of the host galaxies in the optical do- 
main also points to an interaction between an energetic 
(10 54 erg) event and the interstellar medium in galaxies 
(Kulkarni et al. 1998; Wijers 1998). Other evidence in fa- 
vor of such superexplosions is given by Wang (1999), who 
discovered gaseous X-ray emitting remnants of two hyper- 
novae in M101. Some cold HI shells may be the long lasting 
evidence of a short pulse of energy accompanying GRBs, 
as it has been noted by Blinnikov & Postnov (1998) and 
elaborated by Efremov et al. (1998) and Loeb & Perna 
(1998). 

The intention of this paper is to demonstrate the dif- 
ference between HI shells created in connection to GRBs 
and HI shells formed around OB associations. Properties 
of supershells created in connection with a GRB are de- 
termined by a single short-acting source of energy, which 
may arise in any position in a galaxy (and even outside it): 
the energy of 10 51 — 10 54 erg is released in a few seconds 
from a very small volume of the accretion disk around the 
black hole. The activity of young stars in OB associations 
is much less concentrated in space and time: the initial 
volume can have a radius > 10 pc and the energy may be 
delivered during several Myr. The duration of the energy 
supply depends in this case on the history of the star for- 
mation in the OB association and on the IMF, leading to 
a typical spread of 10 - 20 Myr. 

In this paper we compare simulations of expanding 
shells connected with OB associations to those created by 
GRBs. We focus on the star formation around the LMC4 
region in the Large Magellanic Cloud (LMC). Next we dis- 
cuss arguments for and against the origin of multiple su- 
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pershells from single superexplosion; and we consider the 
possibility that GRBs support star formation in galaxies. 

2. Simulations of shells created by a GRB and by 
an OB association 

The expansion of gas layers related to energy inputs from 
GRBs or OB associations can be described as a blast- 
wave propagating into the interstellar medium of the 
host galaxy (Ostriker & McKee 1988; Bisnovatyi-Kogan 
& Silich 1995). Since the radius of the shell is much larger 
than its thickness, the thin layer approximation devel- 
oped by Sedov (1959) can be used. It has been applied in 
one-dimensional models by Chevalier (1974) and in two- 
dimensional models by Tenorio-Tagle & Palous (1987) and 
Mac Low & McCray (1988). These models have been fur- 
ther extended into three dimensions by Palous (1990) and 
Silich ct al. (1996). 

A system of equations of motion, mass, and energy is 
solved. The equation of motion is 



rf_ 

dt 



(mshVsh) — dS sh [{P int - P ext ) + PextVext{Vexp ~ v ext 



+m sh g, 



(1) 



where m s h, v exp , dS s h are the mass, the expansion ve- 
locity and the surface of an element of the shell, Pint and 
P ext are pressures inside and outside the bubble, p ext and 
v ext are the density and velocity of the ambient medium, 
and g is the gravitational acceleration. 
The equation of mass conservation is 



d_ 

It' 



m sh — (v exp - v ext )j_p ext dS sh - m sh T sh . 



(2) 



The first term on the right hand side (rhs), giving the 
increase of mass m s h, is used as long as the normal com- 
ponent of the velocity, (v exp — v ext )±, exceeds the speed 
of sound in the ambient medium. After the expansion be- 
comes subsonic, the mass accumulation stops. The second 
term on the rhs describes the decrease of mass in the shell 
due to mixing with the material inside the bubble; T s h is 
the shell mixing rate. 

The total energy E tot is 



E t ot — E pot + E t h + E ki . 



(3) 



where E pot is the gravitational potential energy in the 
galaxy. E t h and E kin are the thermal and kinetic energies: 



Eth — Eth.int + E t h,sh + E t h ,exti 



-&kin — ^kin.int ~^ ^kin.sh ^kin.ext j 



(4) 



(5) 



where the second parts of the subscripts on the rhs re- 
fer respectively to the medium inside the bubble, in the 
shell, and in the unperturbed environment. In this paper, 



we disregard the thermal energy of the unperturbed ambi- 
ent medium E t h,ext and the kinetic energy of the medium 
inside and outside the bubble Eki n ,int and Eki n ,ext- 

The balance of the internal thermal energy is given as 



—jj_E t h,int — L — Aj„( + T s hE t h,sh, 



(6) 



where L is the energy input rate from a GRB or an OB as- 
sociation and A in t(pint,Ti n t) is the radiative cooling rate 
from the bubble interior, pi nt is the internal volume den- 
sity, Ti n t is the internal temperature. For the mixing rate 
we assume 



sh 



fn 
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where 



fn 



Tint T.sh 
Tint — 



(8) 



T s h is the temperature in the shell. The mixing time t r , 
is estimated as 



tr, 



ar sh 



V, 



sound. sh 



(9) 



where r s h is the radius of the shell and V soun d, s h is the 
sound speed inside the shell: V soun d, s h = & is the 

Boltzmann constant and p is the mean atomic weight of 
particles in the shell, a is a free parameter giving the effec- 
tivity of the mixing process, a lies in the interval (0, oo): 

— a = corresponds to an instantaneous mixing, all the 
mass collected in the shell is immediately mixed into 
the bubble 

— a = oo corresponds to no mixing. 

In our simulations we use a in the range (0.1 - 1), however 
results are almost independent on the particular value, 
and moreover, mixing has influence only on the early evo- 
lution of the supershell. It does not have a substantial 
influence on the later evolution after 1 Myr of expansion. 
The balance equation of the thermal energy in the shell 

is 

-jj.E t h,sh — 2 dt' 1 ( Vexp ~ V ext)\ — h s h — Y s hE t h,sh, (10) 

where the first term on the rhs gives the fraction of the 
kinetic energy of the shell which is converted into the 
shell thermal energy due to the compression of the ambi- 
ent medium. A s h{p s h,T s h) is the cooling rate of the shell, 
where p s h is the volume density in the shell, which is es- 
timated under the assumption that the thickness of the 
shell is a small fraction, (~ 0.1), of the shell radius. T in t 
and T s h are derived from the thermal energies and num- 
ber of particles in the bubble Ni nt and in the shell N s h'- 
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Fig. 1. The shell radius r 3 h as a function of time (left panel) and the expansion velocity of the shell v exp as a function of radius 
(right panel) for shells expanding in a homogeneous medium with p ext = 1.4 cm -3 . Solid lines denote the continuous energy 
input, the dashed ones the abrupt energy input. The thickness of the line gives the input energy, from thin to thick: 10 52 , 10 53 
and 10 54 erg. 



rate from both the bubble and the shell is calculated us- 
ing the cooling functions of Bohringer & Hensler (1989) 
for T € (10 3 - 5 - 10 8 ) K, and Schmutzler & Tscharnuter 
(1993) for T E (10 8 - 10 9 ) K. 

The internal pressure Pint is derived from the internal 
thermal energy E t h,int and the volume of the bubble Vi n t 

„ 2 E t h 

Fint = -z—rj ■ (11) 

-J Vint 

The distribution of the initial energy among the ther- 
mal and kinetic parts corresponds to the solution calcu- 
lated by Sedov (1959), and the division of the initial ther- 
mal energy and mass between the internal volume and 
the shell is quite randomly chosen. However, different ra- 
tios of initial masses and energies do not influence results 
substantially. 

The shell is split into elements in a 3D space, and 
equations of motion ( 1 ) , mas s (p| ) , internal and shell ther- 
mal energy balanc e (|6| ) and (|lfjp, and the equation for the 
internal pressure (|11[) are solved numerically using finite 
timesteps. To keep the accuracy of the integrated quan- 
tities under desired limits, an adaptive step-size control 
scheme is used. 

The expansion in the static, homogeneous medium 
with p ext = 1.4 cm~ 3 is shown in Fig. |l|. We compare 
results of the abrupt and continuous energy input of 10 52 , 
10 53 , and 10 54 erg: 



— In the case of abrupt energy input, the shell rapidly ex- 
pands to the radius of 50 - 150 pc during the first 10 6 
yr. After 2-3 Myr the temperature Ti nt decreases be- 
low 10 6 K. The initial inflation is followed by a slower 
expansion at velocities v exp close to the speed of sound 
in the ambient medium. 

— In the case of continuous energy input, the growth of 
the radius is slower, reaching the size of 50 - 150 pc 
only after 2-3 Myr. There is rapid growth phase as 
in the previous case, however, the expansion velocity 
does not decrease as quickly. Already after 1-2 Myr 
it is higher than in the case of abrupt energy input. 

After ~ 30 Myr, shells reach diameters of 100 - 350 
pc and 300 - 1000 pc and masses of 10 52 - 10 71 and 
10 6 - 9 — 10 8 1 M Q for, respectively, abrupt and continuous 
energy inputs of 10 52 — 10 54 erg. As seen in Fig. [j], shells 
connected to GRBs are smaller than shells connected to 
OB associations for most of the time. At diameters larger 
than 500 pc, GRB shells always expand at velocities below 
10 km s _1 . Thus, large shells with diameters 500 - 2000 
pc expanding at velocities 10 - 40 km s" 1 can not be cre- 
ated by an abrupt energy input connected to GRBs; they 
have to be created by a continuous energy input related 
to massive stars in OB associations. 
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R (kpc) R < k P c > 

Fig. 2. Left panel: a comparison of the observed rotation curve with a model of the LMC composed of dark matter and stars. 
Right panel: the density distribution in a two component disk composed of cold gas - Solid line: (the values of the constant 
density are: 5.0, 1.0, 0.5, 0.1, 0.05, 0.01, 0.001 cm~ 3 ) - supported almost completely by the centrifugal forces, and hot gas - 
Dashed line: (the values of the constant density are: 0.1, 0.01, 0.001, 0.0005 cm -3 ) - balanced almost completely by pressure 
gradients. 



3. Expanding shells in the LMC 

In this section we apply the code described above to ex- 
panding shells in the Large Magellanic Cloud (LMC). For 
the LMC we propose a three component model composed 
of the interstellar medium, stars and dark matter. We as- 
sume that the large-scale shape of the gravitational poten- 
tial is formed by stars and dark matter: the contribution 
of the ISM is neglected. The ISM contributes to local per- 
turbations connected to instabilities and star formation. 

3.1. Stars and dark matter 

The spherically symmetric volume density distribution of 
stars (subscript S) is described by the Hubble-Reynolds 
profile 



Ps(R) 



Ps,o 



(*) ; 



3/2 • 



(12) 



and the less steep distribution of the dark matter (sub- 
script DM) with a function 



Pdm{R) = 



Pdm,o 



1 + (adm) 



(13) 



and Rdm the radial scale length of the dark matter. The 
subscript denotes values at the galactic center. 
The stellar + dark mass within the radius R is 

M S+DM {R) = 4irR 3 s p s ,o x 



4itR 3 dm pdm,o 



arctan(- 



(14) 



The galaxy extends up to the galactocentric distance 



Rb- 

P = Ps + Pdm = for R> R B - 



(15) 



To balance the radial force a circular velocity V rot (R) 
is required: 



V rot (R) = 



^ GMs+DMjR) ^ 1 ' 2 



(16) 



where G is the gravitational constant. 
The potential is given as 



Vrot(R B f 



F(R), 



(17) 



where 



where p is the volume mass density, R is the galactocentric 
distance. Rs is the core radius of the stellar distribution 



F(m = 1 4- 1 r b m dm 
[n > + 2 R DM (M dm + M s ) 
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(RbX 
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DM 
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Rb 



Rdm 

Ms 



R s (M s + M DM ) 



Mi 



R 

Rs 



Rb 
Rs 



fori? < R B 



Rr 



F(R) = -f for R > R B - (18) 

The rotation curve of the LMC given by Kim et al. 
(1998) is reproduced with Rs = 1 kpc, Rdm — 2 kpc, 
Rb = 8 kpc, and the total mass Ms+dm(8 kpc) = 8 x 
10 9 M Q . 60% of this mass is in stars and 40% in dark 
matter (see Fig. ^). 

3.2. Interstellar medium 

The interstellar medium is distributed in a two component 
disk, which is partially supported by the centrifugal force 
and partially balanced by the pressure gradients: 



V G {r) 



e 2 r 



d0 
dr 



1/2 



(19) 



where Via is the plane parallel rotational velocity of the 
gas at a distance r from the rotational axis. A fraction e 2 
of the gravitational force ^ towards the rotational axis 
is balanced by the centrifugal force, the complementary 
fraction (1 — e 2 ) is balanced by the pressure gradients, e is 
a free parameter from the interval 0<e<l:e = l means 
that the gravity is completely balanced by the centrifugal 
force, e = means that the gravity is completely balanced 
by pressure gradients. This model has been described by 
Tomisaka & Ikeuchi (1988) and used by Suchkov et al. 
(1994) and Silich & Tenorio-Tagle (1998) for models of 
large-scale bipolar flows from nuclear starburst in dwarf 
galaxies. 

The distribution of the gas is given as 



PG — Pcold + Phot, 

with 

pi = Pifi exp(const l x), 



(20) 



(21) 



where i stands for cold or hot gaseous components, consti 
is inversely proportional to the square of the sound speed 
cs and to the temperature T of the gaseous component, 
and 



X = <i>{R) - e 2 0(r) - (1 - e 2 )0(O). 



(22) 



The cold and hot gaseous components vary in the degree 
of pressure support: 

— the cold component, where we assume cs = 8 km s _1 
(corresponding to T — 10 4 K), is distributed in a disc 
mainly supported by its rotation: e = 0.9; 

— the hot component, where T = 10 6 K, is more sup- 
ported by pressure gradients: e = 0.5, resulting in an 
almost spherical distribution. 

Densities of both components in a model of the LMC are 
shown in Fig. (§). With p co id,o = 11 cm~ 3 and ph o t,0 = 
0.15 cm~ 3 , the total mass of the cold and hot component 
within 4 kpc is 5.2 x 10 8 M Q , respectively 5.2 x 10 7 M Q . 

3.3. Fragmentation of expanding shells 

Conditions for the instability and fragmentation of ex- 
panding decelerating shells have been expressed by 
Ehlerova et al. (1997) using results of the linear perturba- 
tion theory of transverse motions on a three-dimensional 
shell expanding into a uniform ambient medium, as dis- 
cussed by Elmegreen (1994) and Vishniac (1994). The 
stretching of the perturbed region due to the expansion 
may be compensated by the convergence of the flow due 
to the self-gravity of the shell, if the maximum growth rate 
of a transverse perturbation lo > 0. 




(23) 



where S is the mass column density of the shell, c is the 
speed of sound within the shell. The wavelength of the 
fastest transversal perturbation, A, is given as 



2c^ 



(24) 



io and A can be evaluated in 3D computer simulations us- 
ing the thin shell approximation. In the beginning of an 
expansion, when v exp is large and r s h is small, u) is always 
negative, and the shell is stable. The stretching because 
of the rapid expansion at early stages is much more im- 
portant than the self-gravity of the shell. Later, v exp de- 
creases and r s h increases, so that the first negative term 
in Eq. ( |23| ) becomes less important, and at the same time 
E increases so that the second term in Eq. (^3|), which is 
always positive, becomes more important. At the expan- 
sion time t&, to becomes positive and we start to evaluate 
the fragmentation integral If(t): 



I f (t) - / w(f)df • 



(25) 



At the time t = tf when If(t) = 1, fragments are already 
well developed. We also calculate the wavelength A at the 
times tf, and tf. 
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Fig. 3. R — z cuts of expanding shells at a distance of 2 kpc from the center of the LMC after 10 Myr (left panel) and 20 
Myr (right panel) of evolution. The galactic center is in the direction to the left of the (0, 0) point. The density of the ambient 
medium p ex t{z = 0) from the upper to the bottom row is: 1.4, 10, 20, 30 cm~ 3 . Thin and thick lines correspond to an input 
energy of 10 52 and 10 54 erg respectively, solid and dashed lines mean a continuous or abrupt energy input. 
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Table 1. 3D simulations of expanding shells in the model of the LMC. p ex t is the density of the ambient medium at Ro = 2 
kpc, z — kpc. The value 1.4 cm -3 corresponds to the described model, other values are artificially increased. GRB means an 
abrupt energy input, SN a continuous one. 



3.4- Expanding shells 

We have done computer simulations listed in Table [l| 
The value of p ex t = 1.4 cm -3 corresponds to the density 
at a distance of 2 kpc from the center of the LMC (and 
z = 0). We calculated experiments with higher densities 
of the ambient medium, too. The density could be locally 
increased as a result of sweeping of the mass from other 
places. 

Due to gradients of the ISM density and to the form of 
the gravitational potential, shapes of expanding shells are 
not spherical; they are usually elongated in the direction 
perpendicular to the galaxy plane. This is demonstrated 
in the R — z cuts through expanding shells, which are 
shown after 10 and 20 Myr of expansion in Fig. ||. For 
the continuous energy input only the low energy cases, 
Eq = 10 52 erg, do not blow-out to the galactic halo. For 
Eq > 10 53 erg bubbles blow-out to high z distances. With 
the low (not enhanced) value of the density of the ambient 
medium, the blow-out is radially asymmetrical, as can be 
seen in Fig. |^. This radial asymmetry is the result of gra- 
dients in the distribution of the hot gaseous component, 
which produce the outflow in the direction away from the 
galactic center, as seen at high z-distances in Fig. ||. This 
effect may result in a systematic difference between ob- 
served shapes of supershells at the near and far side of the 
LMC: they should be more round in the north-eastern part 



of the LMC than in the south-western one (for the discus- 
sion of the LMC spatial orientation see Sect. 4). As an 
example the nearly circular supergiant shell SG24 in the 
NE and the elliptical shell SG8 in the SW can be given 
(Kim 1998). The abrupt energy input creates much less 
asymmetrical bubbles, which even for higher energies do 
not blow-out to the galactic halo. These shells stay near 
the galaxy plane even for the highest input energies and 
lowest external densities; they never reach the height of 1 
kpc from the z — plane. This is a more general result, 
connected to the initial rapid expansion and subsequent 
decrease of the shell expansion velocity for the abrupt en- 
ergy input. GRB can lead to a blowout if they explode 
below or above the symmetry plane of the galaxy, or if 
they release high energies (~ 10 54 erg) in thin-disk spiral 
galaxies. 

In Table |l| we give fragmentation times t\> and tf and 
wavelengths of the fastest perturbation, A, given by Eq. 
(p4[). We add simulations with low input energy, Eq = 10 51 
erg, corresponding to the energy of one supernova. For 
this input, energy shells are small and always stable for 
any tested density. It means that an isolated SN explosion 
or GRB of the same energy can not induce gravitational 
fragmentation and trigger star formation. To do so in the 
LMC (and elsewhere) higher input energies are needed. 
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Shells related to continuous energy inputs can frag- 
ment for lower input energies than shells created by abrupt 
ones: in the first case Eq > 10 52 erg is enough for a density 
of the ambient medium in the range p ex t = 1-30 cm~ 3 . 
In the second case Eq > 10 53 erg is needed for p ext > 10 
cm -3 , and E Q > 10 54 for p ext = 1.4 cm~ 3 . 

The time evolution of the fragmentation integral If is 
shown in Fig. [|. In the low density ambient medium shells 
are stable for quite a long time. Except for the continuous 
input of 10 54 erg, neither of the cases reaches i/ before 30 
Myr, which means, that the fragmentation is very slow if 
it happens at all. 

For higher values of the density of the ambient 
medium, the shells fragment. For input energies Eg > 10 
erg, the instability starts earlier for abrupt energy input 
cases, but it proceeds faster in continuous cases than in 
their abrupt counterparts. 

4. Triggered star formation at the edge of LMC 4 
in the Sextant region 

In this section we focus on a region in the LMC near Shap- 
ley's Constellation III, where several OB associations have 
been identified by Lucke & Hodge (1970). A supergiant HI 
shell surrounding the most prominent HI void in the LMC 
can be seen in e. g. Kim et al. (1998) at the same location. 
At the SW edge of this shell, outside the HI void, OB as- 
sociations LH51, LH54, LH60w, LH60e, and LH63 form a 
fraction of the arc. This arc has been named the Sextant 



by Efremov & Elmegreen (1998a), since the OB stars make 
about 1/6 of the complete circle. Its radius is estimated to 
be ^170 pc and the average distance between subsequent 
OB associations in the arc is ~ 32 pc. Near the center of 
the Sextant, two star clusters (HS288 and HS287) have 
been found. The oldest stars in the Sextant are ~ 7 Myr 
old. We assume the age of the two star clusters to be 10 - 
15 Myr. The direct determination of ages of these clusters 
is desirable. 



Here, we would like to discuss once more the scenario 
proposed by Efremov & Elmegreen (1998a), i.e. that ~ 15 
Myr ago the energy has been released from star clusters 
HS288 and HS287 and formed an expanding shell. Later, 
the shell fragmented and triggered the formation of OB 
associations, which during the following 7 Myr moved to 
the presently observed configuration, the Sextant. An al- 
ternative concept is the energy injection from a merging 
event between two compact components of the close bi- 
nary system related to a GRB: the progenitor was proba- 
bly ejected from the nearby massive star cluster NGC 1978 
(Efremov 1998, 1999a; Efremov & Elmegreen, 1998b) as a 
consequence of SN explosions leading to the formation of 
compact components. If this is the case, HS287 and HS288 
have nothing to do with the Sextant. It may explain why 
the two clusters are slightly off center with respect to the 
Sextant arc. Some questions remain unanswered in both 
scenarios: particularly, why we see only a fraction of the 
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circle, and how its orientation relative to the line of nodes 
of the LMC (see below) may be explained. 

We try to reproduce the shape of the Sextant, its ra- 
dius, age and the distances between the OB associations 
using numerical simulations of shells expanding in a model 
of the LMC, taking into account different positions relative 
to the symmetry plane. An alternative scenario is the su- 
pershell expanding in the preexisting HI supercloud, which 
may be a result of the gravitational fragmentation in walls 
of the HI ring surrounding the LMC 4 region. 

The position of the LMC on the sky is specified by the 
line of nodes 8q and the inclination i. We adopt values 
given by Luks & Rohlfs (1992): O = 162°, i = 33°. An 
analysis of proper motions of LMC stars in the Hippar- 
cos catalogue revealed (Kroupa & Bastian 1997), that the 
LMC rotates clockwise. In the combination with the anal- 
ysis of radial velocities by Luks & Rohlfs (1992) it means, 
that the SW half is the more distant part of the LMC and 
the NE half is the nearby one. 

The Sextant itself is almost perpendicular to the line 
of nodes with the radius parallel to it. It is also interest- 
ing, that a continuation of the Sextant arc would cross 
the LMC 4 HI ring nearly perpendicularly. Therefore, its 
radius of 170 pc is almost uninfluenced by projection ef- 
fects, while the deprojected average distance between the 
OB associations in the Sextant is Iob ~ 37 pc. 

Using computer simulations, we search for such initial 
energies and densities of the ambient medium to get un- 
stable shells with radii of 170 pc. With a shell of radius 
170 pc we ask what was the value of the fragmentation 
integral If(t) some 7 Myr ago (7 Myr is the age of the 
oldest stars in the Sextant). We request a value between 
and 1. For small values of If (near 0) the formation 
of stars must have been a very quick process, since they 
formed very shortly after the shell started to be unstable. 
For If near 1 the star formation was somewhat delayed. 
If If — about 7 Myr before reaching the radius 170 pc, 
shells do not fragment early enough and cannot form OB 
stars in the Sextant. If If > 1 7 Myr before reaching 170 
pc, shells also cannot reproduce the Sextant, since stars 
would probably form earlier and the observed arc of OB 
associations would be composed of older stars. 

We ask what is the ratio between the wavelength of 
the fastest transversal perturbation in the shell, A, and 
the shell radius r s h 7 Myr before reaching the radius of 
170 pc. In our opinion, this value did not change since the 
formation of the first stars. Therefore, it should reproduce 
the presently observed ratio of the average deprojected 
distance between the OB associations Iqb to the radius of 
the Sextant: ^ = 0.22. 

The above conditions restrict the possible values of the 
input energy and of the density of the ambient medium: 
— The continuous energy input forms shells of a nearly 

correct size for Eq = 10 53 erg and p ex t — 20 cm~ 3 : 

at t = 4 Myr r s h = 105 pc, // = 0.1, and X/r s h = 

0.18. The radius of 170 pc is reached 11 Myr after the 



beginning of the expansion. Other simulations do not 
meet the conditions so well. 
— The abrupt energy input gives the best solution with 
E = 10 54 erg and p ext = 10 cm" 3 . At t = 3 Myr r sh 
= 130 pc, If = 0.1, and X/r sh = 0.22. The radius 170 
pc is reached at t = 10 Myr. Lower energies or higher 
densities do not fulfill the above conditions. 

We conclude, that for both continuous and abrupt energy 
input we can reproduce the correct size, age and aver- 
age distance between fragments, observed in the Sextant. 
The continuous energy input gives the best solution for 
an input energy of 10 53 erg and density of the ambient 
medium 20 cm~ 3 . In the case of abrupt energy input, a 
higher value of energy (10 54 erg) and a lower density of 
the ambient medium (~ 10 cm -3 ) are needed. 



/ 




-0.2 0.0 0.2 

x (kpc) 



Fig. 5. Projected column density of an expanding shell created 
by the continuous energy input from the expansion center at 
the symmetry plane of the galaxy. North is up and West is 
right, the line shows the direction of the line of nodes, the 100 
pc long arrow points towards the center of the LMC. 

We try to reproduce not only the size of the shell at 
the right time and distances between fragments, but we 
would like to have the correct orientation in the LMC 
and give plausible arguments why the Sextant is only a 
fraction of the total circle. Here, projection effects play a 
role. In Fig. ^ we show the column density in the shell 
with the continuous input energy source of 10 53 erg in the 
symmetry plane of the galaxy at t = 11 Myr, p ex t(z = 0) = 
20 cm~ 3 . The galaxy plane is observed in such a position 
that the line of nodes 9q = 162° and the inclination i = 
33°. The column density is distributed in an elliptical ring. 
Its ellipticity corresponds to projected dimensions of the 
3D egg-shaped shell (see R-z cuts in Fig. ||), not to the 
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Fig. 6. The projected column density of a shell expanding from z = 100 pc above (left panel) or z — —100 pc below (right 
panel) the symmetry plane of the galaxy. North is up and West is right, the line shows the direction of the line of nodes, the 
300 pc long arrow points towards the center of the LMC. 



cosine of the inclination angle as would be the case of the 
projection of the in-plane ring. There are two maxima in 
the column density resulting from projection effects. They 
correspond to those directions where the path along the 
line of sight in the shell is long. Minima appear where the 
line of sight intersects the shell almost perpendicularly. 

To get one-sided structures with one maximum in the 
column density, we calculated models with sources at dif- 
ferent distances above or below the galactic plane. In Fig. 
^ we show two cases with a continuous energy source 
of 10 53 erg placed at z — 100 pc and z = —100 pc, 
Pext(z = 0) = 20 cm~ 3 . The expanding structure is asym- 
metrical because of the asymmetry in the density distribu- 
tion and the gravitational potential. Its most dense part is 
at low \z\ distances near the symmetry plane and its low 
density parts open to high \z\ distances. The projected 
column density shown in Fig. ^| has one maximum and 
it has the shape of an arc. However, the orientation and 
opening of the arc is in both cases (with the source above 
or below the symmetry plane) different from the observed 
position of the Sextant. 

At the edge of LMC 4 the HI density may have been 
increased by sweeping of the mass from the HI void as ob- 
served by Kim et al. (1998). The HI collected there prob- 
ably fragmented and formed giant clouds of increased HI 
density. We assume that the source of energy resided inside 
one of such HI clouds. In Fig. we show results of a simu- 
lation where the shell formed around a continuous energy 
source. The expansion started in a homogeneous spherical 
cloud of density 20 cm -3 and radius 150 pc, where the 
energy source has been placed at an off-center position 




Fig. 7. The projected column density of a shell expanding 
from an off-center position inside a preexisting HI giant cloud. 
North is up and West is right, the line shows the direction of 
the line of nodes, the 100 pc long arrow points towards the 
center of the LMC. 



(about 50 pc from the cloud center). The resulting pro- 
jected column density (FigJ^) does not depend strongly 
on the position of the source inside the parent cloud. The 
energy sources inside a rather large area NW from the 
giant cloud center produce the shape and orientation in 
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the LMC similar to the Sextant. This is shown in Fig. [| 
where the model projected column density is plotted over 
the photographic image of the Sextant area. 




Right Ascension (J3000) 



Fig. 8. A part of the photographic image of the Constellation 
III region including the Sextant arc. The original photograph 
from Boyden Observatory, courtesy of Harlow Shapley, Har- 
vard College Observatory, has been put to computer readable 
form by Knut Olsen, Paul Hodge and Don Brownlee. The pho- 
tographic image is overlayed with isodensity lines giving the 
model column density (the maximum projected column den- 
sity is 4.2 x 10 20 cm" 2 , lines show 90% , 80% ... 10% of this 
value) . 

We conclude that shell forming a fraction of the com- 
plete circle in the correct position can be created by ex- 
pansion which takes place in the NW part of a preexisting 
giant HI cloud at the edge of LMC 4. 

5. Can GRB trigger and support star formation in 
galaxies? 

We have demonstrated that the abrupt energy inputs re- 
lated to GRBs form expanding shells, which can fragment 
and trigger star formation. An inspection of Table [l] shows 
that with energy Eq — 10 54 erg shells fragment in any 
tested density of the ambient medium: p ext = 1 — 30 cm -3 . 
In our opinion, GRB are the likely cause of the supershell 
and the triggering of star formation, particularly when: 

1. there is no stellar cluster inside a supershell or inside 
an arc of young OB stars (providing the normal IMF), 

2. there is no evidence of HVC and/or neighbouring 
galaxies, 

3. there is no hot gas inside the shell, even in the case of 
a young supershell, 



4. there is a system of multiple stellar arcs. 

The hypothesis of Blitz et al. (1999), that HVCs are 
debris after formation of the Local group of galaxies, sug- 
gests large distances to HVCs, which imply for a HVC the 
typical mass of 3 x 1O 8 M of neutral gas and the typical 
size of 25 kpc. In such a case the theory, that the largest 
superbubbles may be created by impacts of HVCs, seems 
to be less plausible. 

Many supershells in M31 (Brinks & Bajaja 1986) and 
M33 (Deul & den Hartog 1990) have a circular form and 
there are no visible OB associations inside them. A num- 
ber of supergiant shells without visible clusters is known 
in isolated galaxies, where the probability of HVC infalls 
is low, for example galaxies NGC 3044 and NGC 4631 (see 
references in Efremov et al. 1998 and Loeb & Perna 1998). 

Recently, a special photometric search for clusters, 
which could be the energy source of the supershells in 
the irregular dwarf galaxy Holmberg II, was carried out 
by Rhode et al. (1999). They found only 4 out of 50 su- 
pershells with clusters of suitable ages, which could have 
contained enough OB stars to form the observed super- 
shells, assuming the normal IMF. The best examples of 
supershells with no clusters are the most energetic ones 
found in the low density environment at the periphery of 
Ho II, where no star formation occurs and where the pres- 
ence of massive clusters is almost excluded. Another result 
given by Rhode et al. (1999) is that nearly no HVCs exist 
around Ho II, which may eventually be responsible for the 
observed supershells. Hypernovae are also not the likely 
energy source, since they reside near the star forming re- 
gions (Paczynski, 1998), and there is no star formation 
near the supershells at the periphery of Ho II. The en- 
ergy source could be a merging event of a compact binary 
system. Supernova explosions, leading to the formation of 
compact components, may give the system the high space 
velocity of >100 km s" 1 (Blaauw 1961). Then, the GRB 
occurs far from the formation site as demonstrated e.g. by 
Lipunov et al. (1997). 

The Sextant arc is one in the system of three stel- 
lar arcs in the LMC 4 region (Hodge 1967, Efremov & 
Elmegreen 1998a). If the GRB theory of the origin is cor- 
rect, then the progenitor of the Sextant arc (and progen- 
itors of the other arcs as well) probably belonged to the 
massive 2 Gyr old cluster NGC 1978 located 0.5 - 1.0 
kpc from its center. NGC 1978 has an elongated shape 
which may be the result of merging of two clusters, or of 
disk-shocking. Such events are able to eject many stars 
(including compact binaries) from the cluster. Actually, 
there is concentration of X-ray binaries in the discussed 
region. Also, one probable relict of the merging event of 
the compact binary system, the Soft Gamma Repeater 
SGR 0526-66 in SNR N49, is at a distance of 18' from the 
cluster NGC 1978 (Efremov 1998, 1999a, b). 

At least three systems of multiple stellar arcs are 
known in the LMC, NGC 6946 and M 83 (Hodge 1967; 
Efremov 1999a). The opening angle in most of these arcs 
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is few dozen degrees. The strictly circular form of these 
arcs implies an origin from the central source. Regular 
distances between clusters along an arc suggest their for- 
mation by the gravitational instability in the swept-up gas 
ring (Sextant and two arcs in M 83). 

Many questions related to the origin of supershells and 
triggered star formation in connection to GRBs remain 
open: 

— Can isolated galaxies with active star formation be 
those in which there were recently many GRBs? 

— Can GRBs support and revive star formation activity 
together with OB stars and SNe? 

— Can GRBs be the initial source to the star formation 
and turbulent motions in gas disks of galaxies? 

The last item (Hodge 1998, private communication) is 
quite attractive: compact binaries with neutron stars may 
form in the process of dynamical evolution of massive and 
old star clusters. The two compact components merge and 
produce GRBs long after they escaped from the cluster. 
Those, which explode within the gas disk, might trigger 
the turbulent motions and star formation there. 

Considering the large uncertainties in the theory of 
GRBs, properties of supershells and stellar arcs suggested 
to be created by GRBs may give essential constraints to 
the GRB theory. The central angle of stellar arcs, for ex- 
ample, may be connected somehow with the beam angle 
of the GRB - especially if they cannot be explained as pro- 
jected shells created by sources outside the galaxy plane 
or by sources inside giant interstellar clouds. 

Our simulations show one possible obstacle to GRBs 
playing an active role in star formation events in galaxies. 
The energy needed to create gravitationally unstable shells 
by an abrupt explosion must be very high (E n > 10 53 
erg; for low ISM densities even E > 10 54 erg). There 
were indeed GRBs with energy 10 53 and 10 54 erg, such 
as GRB 971214 and GRB 990123 (Kulkarni et al. 1999), 
but the beaming of the relativistic jets would imply that 
the emitted energy would be proportionally smaller. So, if 
the energy loss during the time interval between the burst 
itself and the time when the expanding shock forms, is 
large, or if the energy released by the burst is small, then 
the effect of GRBs on the host galaxy is minor and it can 
not induce star formation. 

6. Conclusions 

We discussed properties of HI shells connected either to 
GRBs or OB associations: in the first case, there is an 
initial inflation or rapid growth of the shell followed by a 
slow expansion. In the second case the initial expansion is 
much slower, however, the expansion velocity decelerates 
to low values longer. Supershells of diameters larger than 
500 pc are almost static if related to GRBs. If a supcrshell 
larger than 500 pc expands with a velocity higher than ~ 
10 km s _1 , it is undoubtedly related to an OB association 
(= continuous energy input). 



Using computer simulations of expanding shells in 
galaxies, we showed that both the abrupt energy input 
connected to a GRB or a continuous energy input con- 
nected to an OB association can trigger star formation. 
Shells related to GRBs need the higher energy to frag- 
ment into protoclusters than shells connected to OB as- 
sociations, where the lower energy is sufficient. The star 
formation triggered by GRBs can initiate or complement 
the star formation triggered by OB associations. 

As an example, we discussed the Sextant region in the 
vicinity of LMC 4. This arc of OB associations can be 
interpreted as triggered star formation within the HI shell. 
We found that both abrupt and continuous energy inputs 
can correctly reproduce the radius, the age of the arc and 
the average distance between OB associations. 

Partial shells formed by an energy source above or 
below the galactic plane are projected to the tangential 
plane perpendicular to the line of sight as arcs with apexes 
turned to or from the line of nodes of the LMC. This is 
not the case with the arcs in the LMC 4 region. To repro- 
duce the orientation and position of the Sextant we have 
to assume, that the expansion took place inside a giant HI 
supercloud, which had been previously formed at the rim 
of LMC 4. 

Clusters HS287 and HS288 are not in the very center of 
the Sextant arc. This may be explained as a consequence 
of the asymmetry of the shell: the projected structure (or 
its densest parts) has a center at a position different from 
the energy source. 

A comparatively large number of SNe is needed to fit 
the observed arc. A cluster producing 10 - 100 SN requires, 
with the normal IMF, many low-mass stars in the same 
cluster. Unfortunately, there is no photometry for these 
clusters available to check their mass function. Because of 
the proximity of the massive star cluster NGC 1978, there 
is a high probability that both compact binaries and hy- 
pernovae exist in this region (e.g. Chevalier & Li 1999). 
The occurrence of very massive clusters near the arc sys- 
tems in LMC and NGC 6946 shows that their progenitors 
may be stellar objects. 

In any case, the presence of two or three more arcs 
composed of OB stars and clusters in the same region near 
LMC 4 (Efremov & Elmcgreen 1998a), and the absence of 
similar arcs elsewhere in the LMC, can be better explained 
by a number of progenitors escaping from a single massive 
cluster, presumably NGC 1978. The recent discovery of 
the most massive young cluster in the galaxy NGC 6946, 
which is located near the only multiple arc system in this 
galaxy (Hodge 1967; Larsen & Richtler 1999), supports the 
idea of the connection between multiple, not concentric 
arcs with progenitors escaping from one common massive 
cluster. 
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